Violation of action-reaction and self-forces induced by nonequilibrium fluctuations 
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We show that the extension of Casimir-hke forces to fluctuating fluids driven out of equilibrium can exhibit 
two interrelated phenomena forbidden at equilibrium: self-forces can be induced on single asymmetric objects 
and the action-reaction principle between two objects can be violated. These effects originate in asymmetric 
restrictions imposed by the objects' boundaries on the fluid's fluctuations. They are not ruled out by the second 
law of thermodynamics since the fluid is in a nonequilibrium state. Considering a simple reaction-diffusion 
model for the fluid, we explicitely calculate the self-force induced on a deformed circle. We also show that the 
action-reaction principle does not apply for the internal Casimir forces exerting between a circle and a plate. 
Their sum, instead of vanishing, provides the self -force on the circle-plate assembly. 
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In his pioneering work, Casimir [1] showed that two metal- 
lic plates in the electromagnetic field of vacuum attract one 
another due to the restriction they impose on the quantum fluc- 
tuations of the field. Fluctuation-induced forces exerting be- 
tween macroscopic objects have since been exhibited in a vast 
variety of other systems at equilibrium [2], such as critical flu- 
ids and crystal liquids in the nematic phase, in which thermal 
fluctuations of long range can develop [3, 4], These forces are 
thought to play an important role in the stability of equilib- 
rium phases of mesoscopic particles embedded into complex 
fluids. While the theoretical achievements in this field are nu- 
merous, a direct measurement of the Casimir force in critical 
fluids has only been obtained recently [5]. 

Matter in nonequilibrium steady states also develops fluc- 
tuations that are generically of large amplitude and have long 
correlation lengths [6, 7]. By analogy to the equilibrium sit- 
uation, one expects these fluctuations to induce similar forces 
that may be responsible for the aggregation and/or segregation 
mechanisms observed in fluids driven out of equilibrium [8- 
10]. However, the calculation of these forces cannot rely on 
the derivation of an (equilibrium) thermodynamic potential. It 
is only recently that these forces have been obtained between 
two planar objects immersed into nonequilibrium driven sys- 
tems [11], granular fluids [12], or reaction-diffusion systems 
that violate the detailed balance [13]. 

Whether at or out of equilibrium, accurate experimental 
measurements of fluctuation-induced forces need to go be- 
yond the idealized geometry of infinitely-long plates, predom- 
inant in theory for its simplicity. Although a long-studied 
topic, the proper account of the geometry dependence of 
Casimir forces is a notoriously difficult problem when deal- 
ing with nontrivial geometries [14]. To date, the most widely 
used technique by experimentalists [5] relies on the so-called 
Derjaguin construction [15] (proximity force approximation) 
which in essence integrates the two-plate expression of the 
force along the curved surfaces. 

In this Letter, we show that when restricted by nonplanar 
objects, nonequilibrium fluctuations can lead to additional ef- 
fects not possible in equilibrium systems. Namely, nonvanish- 
ing forces can be induced on single asymmetric obstacles, and 



the action-reaction principle between two intruders can be vi- 
olated. These phenomena can have significant consequences 
in experiments, as they lead to directed motion and uneven- 
ness in the measures of the forces between two objects. Fur- 
thermore, an unbalance of action-reaction would impeed the 
use of the Derjaguin approximation where it would normally 
be valid at equilibrium. We found no mention of these facts in 
the literature. 

Since nonequilibrium systems are thermodynamically 
open, it is not entirely unexpected that self-forces can appear. 
In fact, provided that both the microscopic time-reversibility 
and space rotation-invariance symmetries are broken, such 
forces have been implicitely suggested by the occurence of 
sustained motions in other nonequilibrium contexts, such as 
in ratchets [16], Brownian motors [16, 17], molecular mo- 
tors [18], or the adiabatic piston [19]. In these systems, the 
space asymmetry usually lies in an external temperature gra- 
dient or in an anisotropic field exerting on the object. More 
recently, however, the directed motion of an asymmetric ob- 
ject immersed into vibrated granular matter has been exhib- 
ited [20]. The direct calculation of self- forces that we present 
here allows for a better understanding of the different effects 
at play in such motions. It also makes possible the evalua- 
tion of additional stresses exerting on asymmetrical structures 
in micro-devices, and could be used in tailoring mechanisms 
for the self-assembly of ordered structures. The violation of 
action-reaction between two intruders directly results from 
the presence of self-forces. It does not seem to be system- 
atic, however, even in asymmetric setups [21]. Note that it 
prevents the two-body forces to derive from an effective po- 
tential, in contrast to equilibrium cases. Let us add that a vio- 
lation of Newton's third law has also been noted in depletion 
forces between identical spheres in a flowing fluid [22]. 

For illustration purposes, we exhibit here these effects with 
the rather simplified nonequilibrium fluid that has been used 
in [13] in a planar geometry. The applicability of the model 
to nontrivial geometries can be greatly facilitated by devising 
an adequate Green function formalism. We then show that the 
self-force exerting on a deformed circle is indeed nonzero at 
second order in the radius perturbation when dipolar defor- 



mations are considered. We also calculate the internal forces 
between a circle and a plate in an asymptotic regime where 
the circle's radius is small in comparison to the correlation 
length while its separation to the plate is large. In this situa- 
tion, action-reaction is not satisfied as the circle-plate assem- 
bly experiences a net self-force. 

The fluctuating media is described by a reaction-diffusion 
fluid, whose nonequilibrium steady state is achieved by vio- 
lating the detailed balance [23]. The local density p{r,t) of 
the fluid fluctuates around a homogeneous reference density 
po = (p(r,f)), where (•) is a stochastic average. The density 
deviation <I>(r,f) = p(v,t) — po is assumed to satisfy, at the 
mesoscopic scale, the stochastic reaction-diffusion equation 



~dt 



= (DV2-7)4> + <^, 



(1) 



where D is a diffusion constant, 7 the reaction rate that drives 
the system to local equilibrium, and <^ (r,f ) is a random white 
noise of correlation intensity Y that takes into account the fluc- 
tuations on the reaction rates. Eq. (1) primarily describes 
density fluctuations in a fluid with two reacting and diffus- 
ing chemical components (see [23]), but other nonequilibrium 
systems in their steady state are described by this model. The 
steady state fluctuations are characterized by the bulk correla- 
tion length K^' = (7/D)^''^ that can be chosen as the meso- 
scopic scale. 

Static objects immersed in the fluid prevent any flow of mat- 
ter across their surface. Eq. (1) is thus supplemented by the 
non-flux condition n • VO = at the objects' surface, where 
n(r) is a unit normal vector pointing outward from the fluid's 
domain. In a steady state, one expects the pressure p of the 
fluid to be related to the density by a local equation of state 
P = /'(p(i',f))- This relation is experimentally measured in 
a number of cases of interest, like in driven granular media 
[24], for example. Here we only assume that it is expandable 
around the reference density po and that density fluctuations 
stay small. The average pressure is thus modified by the fluc- 
tuations according to {p) = pQ + ^ (*^^) where pQ — p{po) 
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and Pq ~ 3-7 (Po)- The total force F5 exerted by the fluid on 
an immersed object S results from summing this local pressure 
on every element da of its surface. Since po is a homogeneous 
pressure, it does not induce a force and one has 



F, = ^^dcTn(4>2). 



(2) 



Starting from a dynamical model like (1), Casimir forces 
can be obtained by using Green functions [4]. After an initial 
transient of characteristic time of order 0(7^'), the stationary 
solution of (1) is 

^,i{r,t) ^ JAt' j Ar' G{ry ,Dt-Dt')^{r' ,t'), (3) 

where i2 is the domain occupied by the fluid. 

It can be established that the temporal Fourier transform of 
the Green function G(r,r',(B) ~ /dTe"*'^G(r,r',T) with T — 



Dt, is, up to a factor, the static structure factor of the fluid that 
enters into the force (2) when evaluated at o = 0: 

(4>,(r,f)4>,(r',f)) = ^|g|^dr"G(r,r",a))G(r",r',-a)) 



= — G(r,r',a)-0). 



(4) 



The first equality in (4) follows directly from (3) and the 
convolution theorem. The second results from the differential 
equation satisfied by G(r,r', (o) (see details in [25]): 

(- V2 + K-2 - /a))G(r,r', «) = 5(r - r'), (5) 

n(r)-VG(r,r',co)|rGS = Vr' e ii,Va). (6) 

In view of (4) one can omit any further reference to © and cal- 
culate G(r,r') solution of (5)-(6) in which (O is set to right 
away. This is an appreciated simplification; F5 only depends 
linearly on the static Green function. 

To deal with the difficulties brought about by nonplanar ob- 
jects in Casimir forces, a natural way is to use multiple scat- 
tering techniques [14]. If Go is the free space (unconstrained 
fluid) Green function, (5)-(6) are also equivalent to [25] 

G(r,r') = Go(r-r')-ydcJi G(r,ri) ni-ViGo(ri-r'), (7) 

which we abbreviate as G = Go + G* VGo. The recursive 
iteration of this integral equation expands G as a series of 
multiple scatterings of Go on the surface S. In three di- 
mensions, Go(r) = exp{— K'|r|}/47r|r| and in two dimensions, 
Go(r) = Ko(K'|r|)/27r, where Ko is the modified Bessel func- 
tion of order 0. As it is obvious from (7) and the above ex- 
pressions for Go, the problem of calculating the force (2) as 
it stands is ill-defined: short-range divergences appear. They 
are due to the inaccuracy of the continuous model on the mi- 
croscopic scale [13]. A "bulk" divergence occurs when evalu- 
ating G at a same point r, as well as a "wall" divergence once 
this point is approached to a surface. The first divergence is 
trivial to remove: it is independent of the immersed objects 
and thus consists in a homogeneous (although infinite) pres- 
sure unable to produce a force. The wall divergence plays 
an important role: it originates in the boundary condition im- 
posed on G(r,r') by the immersed objects and it is integrated 
all along their surface in calculating the force. The issue is 
then to understand how this integrated divergence compen- 
sates itself between different sides of the objects to yield a 
finite result. This compensation does not occur for any shape. 
To illustrate this, we consider the Green function Gp of a 
fluid restricted to a half space by a plate. The condition (6) 
on the plate can be replaced by the addition of an "image" 
source 5(r — r'*) on the r.h.s. of (5), where r'* is the point 
symmetric to r' w.rt. the plate. The solution hence reads 
G/>(r,r') = Go(r - r') + Go(r - r'*). It is clear that evalu- 
ating {Gp — Go)(r,r) (having subtracted the bulk divergence 
Go(r — r') as r' -^ r) at the plate's surface produces a diver- 
gent collapse of r and its image r* . For a smooth surface, 
this divergence is only slightly modified by the curvature and 



one expects the total force to be finite. By contrast, objects 
with sharp corners produce several images of r (for instance, 
three in a right corner in 2D). They generate additional diver- 
gences in the edges that are not likely to be compensated (un- 
less a symmetric corner exists on the other side of the object). 
We restrict ourselves to objects with radii of curvature large 
enough to avoid such complication. The proper mathematical 
way of removing the divergences is to introduce short-range 
(e.g., "hard-core") cutoffs that are removed at the end of the 
calculation, similar to classical Coulombic systems between 
opposite charges and at metallic walls. Using (4), the regular- 
ized force Fj (2) on the object S is thus given by 

Fs - \imFoK /dcjn(r) [G-Gq] (r-en(r),r-en(r)), (8) 

£^0 .Js 

where Fq = PqF / {4Dk) has the dimension of a force. We 
now apply this general framework to calculate the fluctuation- 
induced force (8) on two distinct systems embedded into the 
fluid. For simplicity, we limit here the fluid to two dimensions, 
but the conclusions also apply to three-dimensional cases. 

Deformed circle. Since the force (8) on a single circle 
must vanish by symmetry, we deform its radius R accord- 
ing to R{9) = R + ris{0) (in polar coordinates), and assume 
7] <C JC^^ R. The general solution of (5) for a finite object is 



G(r,r')=Go(r-r') + y: 






In 



-amn^miKp)K„{Kp'), (9) 



where (p,0) are polar coordinates for r, {p',9') for r', and 
Km is the modified Bessel function of order m. The coeffi- 
cients Omn satisfy «„,„ — a„m = al^ _„ to ensure that the Green 
function is real and symmetric under the interchange of r and 
r'. They still need to be determined from the boundary condi- 
tion at p = R{9). Substituting (9) in (6), one can obtain them 
perturbatively in 77 in terms of the Fourier coefficients ofs{9), 
s„ = i2n)-'J^''d9e-"^s{9). 

In the force (8), it must be noted that the dependence on 
R{9) is double: explicit in da, n, and G (via flm„), and im- 
plicit since the Green function is evaluated on the boundary of 
the deformed circle. The whole expression is expanded in 77 
and it is verified that the zeroth order contribution, which cor- 
responds to the force on the undeformed circle, vanishes. The 
contribution linear in the perturbation also vanishes. Indeed, 
by linearity, each Fourier mode of s can be analyzed sepa- 
rately; the mode n =0 merely corresponds to a change in the 
radius of the circle; the dipolar modes n ~±l are equivalent 
to a small displacement of the unperturbed circle; all remain- 
ing modes |«| > 2 correspond to symmetric perturbations, not 
having any preferred direction. The first contribution to the 
force thus comes from the second order in 77 and has the form 



F = FQ{Kri)^Y,f'"^>"S- 



X—m-) 



(10) 



where the two components of F are expressed in the r.h.s. 
as a complex number. This particular combination of the 
(complex) Fourier coefficients i„ is necessary to ensure that 



the force correctly transforms as a vector, or, equivalently in 
Fourier modes, as a dipole. The real coefficients /„, are non- 
trivial series of Bessel functions evaluated at /? or /? + e, and 
the limit e ^ must be taken after summing them. 

For simplicity, we choose a specific deformation that leads 
to a nonvanishing force. It can be checked that fo = f-i =0, 
as expected because the joint deformation given by ^o and ^i 
produces again a circle that is merely expanded and translated: 
its self-force vanishes. The first nontrivial cases are given by 
the coupling of the n = 1 and n = — 2 modes of s{9). Hence 
considering s{9) — 2sicos{9) +2s2cos(20), we find 



^FQSiS2iKri)^H{KR)x, 



(11) 



where // is a dimensionless function whose numeric com- 
putation is accurately compatible with 2/(kR) in the range 
0.1 <kR< 10. 

A nonvanishing self-force is therefore produced on a de- 
formed circle at order 7j^ for the simple deformation consid- 
ered here, made of a dipolar and quadrupolar combination. 

Circle-plate system. In systems at equilibrium, 
fluctuation-induced forces can only appear between two 
or more objects. Because the global force on the total 
system must vanish, such two-body forces always satisfy the 
action-reaction principle. The picture is different in systems 
driven out of equilibrium. 

Indeed, consider two objects S and S' immersed in the 
fluid. The total force F5 on S can be separated into a self 
contribution F^, that may already be present in the absence 
of 5', and a contribution F^^^/ = F5 — F§ due to the addi- 
tional asymmetry provoked by the presence of S' . Denoting 
by G^y and G^ the Green functions associated to the two- 
object and single-object S systems, respectively, one has from 
(8) ¥s^s' ^FQKJsdan{Gss' - G^) and 



-s-^s' 



-s'-^s - 



^55' 



'/ — r c — J 



(12) 



where F^y is the global force exerting on the assembly 5 U 5' 
considered as a whole. Note that the quantity G^y — G^ en- 
tering in F^^y is well-defined: the bulk divergence and the 
wall divergence on S, present in both Gss' and G^, compen- 
sate in the subtraction. We take as definition of the action- 
reaction principle for the internal (two-body) forces of such 
system the vanishing of the r.h.s. of (12). Since fluctuation- 
induced forces are not additive, this vanishing will not happen 
in general in the presence of self-forces, except from obvious 
symmetry reasons. 

As an example, we take an assembly made of a circle C of 
radius R and a thin and long plate f in a 2D fluid. Their sepa- 
ration at the closest point is d (see Fig. 1). The plate is taken 
much longer than K"^' to avoid boundary effects. Since either 
object is symmetric, F[1 = F^ = 0. We calculate both terms 
in the l.h.s. of (12) in the regime R <^ K^^ <t^ d to show that 
the total force on the assembly, Fcp, is nonzero; equivalently, 
action-reaction is not satisfied in this situation. 

In the multiple-scattering expansion of Gcp(r,r)\rec (7), 
the free-space correlation Go is scattered on both surfaces C 



e 



FIG. 1 : A circle C of radius R and a thin, long, plate P are immersed 
in the fluctuating fluid. Their separation at closest point is d. 



and P. It is clear that when the separation d is much larger 
than the correlation length fc^', the dominant terms in this ex- 
pansion contain the least number of propagations between C 
and P. Those containing none rebuild the series of scatter- 
ings of Go on C that gives the Green function of the single 
circle, G^. All other terms contain scatterings on P. Since 
the initial and final points are both r EC, making a scattering 
on P necessarily involves at least two interspace propagations, 
one from C to P, and one to come back to C. The next dom- 
inant terms therefore contain exactly two of the propagators 
Go evaluated at points separated by at least d. However, be- 
fore and after these interspace propagations, any number of 
scatterings from C to C or from P to P can be done. Summing 
them up into the quantities G^ and G^, one eventually has 



(Gcf-G«)(r,r)|rec''- 



'G?.*VG^- 



-Gc*yG°P*yGl. 



Explicit expressions for GJI and G^ are easy to obtain. The 
function G^ is straightforward to calculate from (9) with the 
result fl„„ ^ -{V„{kR) /K'„{KR))8n-m, where I^ and K;„ are 
the derivatives of the modified Bessel functions of order m. 

We shall here only state the result of an asymptotic anal- 
ysis of Fc^p and F/>^c, based on small- K"/? expansions and 
steepest-descent values of integrals asKd ~>°° (explicit calcu- 
lations, including in other regimes, will be developed in [25]). 
On the X axis, the forces exerting on the circle and the plate in 
the regime R <^ K^^ ^d are found to be 



Fc. 



-Fo 



^/n{KRfe 



2^-2Kd 



Kd 



Fp. 



-^Fc^P. (13) 



Clearly, action-reaction is not satisfied. Furthermore, the 
circle-plate assembly experiences a nonzero global self-force 

along X of the same order: F^p ^ \Fo v'^t'^ )f 

In conclusion, we have shown that in nontrivial geometries 
the extension of Casimir-like forces to fluctuating fluids driven 
out of equilibrium must take into account two interrelated phe- 
nomena forbidden at equilibrium: the possibility that a self- 
force may be induced on a single asymmetric rigid body, and 
that the action-reaction principle for the forces between two 
objects may be strongly violated. The latter fact impeeds that 
an effective interaction potential holds in nonequilibrium. Its 



occurence would prevent the use of the Derjaguin approxima- 
tion. As the magnitude of this violation can be of the same 
order as the internal forces, special care should be taken when 
obtaining these forces in experiments or simulations. 

The complexity of dealing with nonplanar objects has been 
overcome here by considering a very simple model for the 
nonequilibrium fluid and by devising a Green function and 
multiple scattering technique. Clearly, to allow quantitative 
comparison with real fluids (such as colloidal solutions, dusty 
plasmas, etc.), one would need to refine the initial model. 

The presence of self-forces leads to directed motion if the 
objects are let free to move as in the case of ratchets. Self- 
forces can lead to arrangements of composites of intruders that 
could be tailored to produce microdevices by self-assembling. 
Their dynamical properties, however, need a more thorough 
analysis, for their motion will affect the fluid's fluctuations 
and a self dynamical interaction could appear. 
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